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We present a detailed study of a linear sigma model containing one chiral nonet transforming 
under U(1)a as a quark-antiquark composite and another chiral nonet transforming as a diquark- 
anti diquark composite (or, equivalently from a symmetry point of view, as a two meson molecule). 
The model provides an intuitive explanation of a current puzzle in low energy QCD: Recent work has 
suggested the existence of a lighter than 1 GeV nonet of scalar mesons which behave like four quark 
composites. On the other hand, the validity of a spontaneously broken chiral symmetric description 
would suggest that these states be chiral partners of the light pseudoscalar mesons, which are two 
' quark composites. The model solves the problem by starting with the two chiral nonets mentioned 

, and allowing them to mix with each other. The input of physical masses in the SU(3) invariant limit 

■ for two scalar octets and an "excited" pion octet results in a mixing pattern wherein the light scalars 

have a large four quark content while the light pseudoscalars have a large two quark content. One 
light isosinglet scalar is exceptionally light. In addition, the pion pion scattering is also studied and 
the current algebra theorem is verified for massless pions which contain some four quark admixture. 



00 

o 



O 

o 



■ 

X 



PACS numbers: 13.75.Lb, ll.15.Pg, ll.80.Et, 12.39.Fe 



Dh; I. INTRODUCTION 

D , The topic of anomalously light scalar mesons in QCD has become a subject of increasing interest in the last fifteen 
years or so [H-fl^- Of course light scalars, especially the light isoscalar called first sigma and now /o(600), have been 
discussed for at least three times as long, although without general agreement on their actual existence. The difficulty 
—i people had previously in accepting the light scalars was largely due to the great success of the simple quark model, 
^ in which the lightest scalars are expected to be p-wave quark- antiquark composite states and hence to be in the 1 
CO to 1.5 GeV range, like the other p-wave states. It seems that physicists now believe more in their existence because 
there have been an increasing number of investigations, using a variety of techniques and models, which suggest that 
they do exist. Common features in many of these approaches have been the use of unitarity (which no one denies) 
and some input at low energy from chiral dynamics (which is also considered reasonable) . 

Of course, the strongly interacting gauge theory QCD has not been "solved" and any possible new features in 
the low energy region where the effective coupling constant is especially strong raise the hope of improving one's 
understanding of this basic theory. Perhaps the most fascinating possibility is that the very light scalars contain two 
quarks and two antiquarks. Variants based on a diquark- anti diquark picture [4^1 or a meson-meson "molecule" 
picture ^if have been discussed. 

A lot of attention has been given to the question of a possible nonet grouping for the light (less than 1 GeV) scalars. 
■ The candidates are the already mentioned /o(600), the Kappa(800-900) [not conclusively estabhshed according to [l|n, 
Ci I the established ao(980) and the estabhshed /o(980). It has been pointed out [See for examples [l^l, [l3| and p^ ] 
that a characteristic signature of a four quark content would be an inverted mass ordering, with an almost degenerate 
/ = 0, / = 1 pair being the heaviest rather than the lightest states when the light quark masses are " turned on." 
This seems to be the case. 

Associating the four quark states with the lightest scalars naturally raises the question of where are the p-wave quark, 
antiquark scalars. The candidates for the non-zero isospin states are the established ao(1450) and the established 
K^{1430). For the / = states the established candidates are the /o(1370), /o(1500) and the /o(1710), one of which 
may be a glueball. There is a slight puzzle since the non strange ao(1450) with a listed mass of 1474 MeV is heavier 
than the strange Kq{1A30) with a listed mass of 1414 Mev. In addition some branching ratios are not well predicted 
by SU(3) invariance. A possible way to overcome this problem [[29|, |33|] is to allow mixing between the lighter 4 
quark and heavier 2 quark scalar nonets. This feature is incorporated as a basic part of the present paper. 
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While it is rather difRcult to treat low energy QCD dynamically, great success has been obtained at low energies 
using the underlying chiral symmetry of QCD. We also incorporate this as an aid in getting more information about 
the system. This feature will be implemented by using linear rather than the more usual nonlinear representations 
for the pseudoscalar and scalar fields. If both four quark and two quark scalars are present, this means that four 
quark pseudoscalars should also appear in the model. Experimental candidates for the non-zero isospin, higher mass 
pseudoscalars are the 7r(1300) and the two not yet conclusively established strange states ii'(1460) and K{1830). The 
candidates for the higher mass isoscalar pseudoscalars are the 77(1295), ?7(1405), ?7(1475) and the not conclusively 
established 77 (1760). It is possible that one or more of these experimental candidates also contain glueball and radial 
excitation admixtures. At first glance it might seem puzzling that the picture seems to be: light mass two quark and 
heavy mass four quark states for the pseudoscalars at the same time as light mass four quark and heavy mass two 
quark states for their "chiral partners," the scalars. We shall make no initial assumption on this matter but let the 
experimental particle spectrum together with the mixing inherent in the model tell us the answer. 

At a technical level it is amusing to note that the U(1)a transformation properties distinguish the two quark from 
the four quark fields. Since it is known that the U(1)a symmetry is badly broken in QCD, this means that we have to 
model the breaking in some detail. For this purpose we will use an extra term in addition to the usual one. We adopt a 
counting scheme for selecting the most important terms, out of the many possible ones. We assign a number N equal 
to the number of underlying quark plus antiquark lines associated with each effective term. Then it seems reasonable 
to pick up the terms with smallest N values. On this basis the extra term for saturating the U(1)a anomaly has the 
same justification as the conventional one. 

Clearly, with so many scalar and pseudoscalar fields present, the model is fairly complicated to analyze. At the same 
time it is widely believed that massless (ie zero mass light quarks) QCD is an excellent qualitative approximation. 
Except for the pseudoscalar Nambu Goldstone bosons of the theory, the masses of the physical particles made from 
light quarks are largely due to the spontaneous breakdown of chiral symmetry. We will employ this limit of the theory 
in the present paper and note that it very much simplifies the analysis. Especially, the characteristic mixing matrix 
pattern of the two quark and four quark states becomes very clear. The puzzle of opposite two quark vs four quark 
structures of the scalars and pseudoscalars seems to be neatly solved by the mixing mechanism. 

Even though the nonlinear chiral model is more convenient for systematically studying the loop corrections at 
very low energies, the linear sigma models have a long history of elucidating key features of the strong interactions. 
Roughly speaking the use of the nonlinear model amounts to integrating out the scalars (although it is technically 
somewhat more general than that). Certainly for learning about the scalars themselves it is rather convenient to have 
them present in the Lagrangian to begin with. The most famous example of the linear model is of course the Higgs 
potential of the standard model. One of the classical triumphs of the nonlinear model is the derivation of the "current 
algebra" formula for low energy pion scattering. This can be obtained, though in a more complicated way, also in the 
linear model. We verify this in detail in the present paper. One might wonder, since the pion in the present model has 
a small (but non negligible) four quark content, whether the current algebra result strictly does hold in the present 
model. Our result shows that it does hold for the zero pion mass case we are considering here. 

The two chiral nonet model was introduced in [sij as a convenient way to study the possibility of mixing between 
quark- antiquark (qq) spin zero mesons and two quark- two antiquark (qqqq) spin zero mesons. Altogether there are 
two pseudoscalar and two scalar nonets contained in the model. It was found that, in the zero quark mass limit 
with just a few explicit chiral invariant terms contained, there was a possibility of a situation in which the lightest 
pseudoscalars could have zero mass (i.e. be Nambu- Goldstone bosons) and be primarily qq type while the next heaviest 
mesons could be scalars, primarily of qqqq type. Furthermore, the next heaviest mesons could be pseudoscalars of 
mainly four quark type while the heaviest could be scalars, mainly of two quark type. A treatment [43| of the model 
with similar chiral invariant terms and several different quark mass terms also found that light scalars with relatively 
large admixtures of qqqq type states are favored. Actually, the model can be rather complicated since there are twenty 
one renormalizable chiral invariant terms which can be made as well as a similar number of renormalizable quark 
mass type terms which transform as the (3, 3*) -I- (3*, 3) representation of chiral SU(3)lX SU(3)r. In the present 
authors studied the more general version of the model in which all possible chiral invariant, even non renormalizable, 
terms were included together with the single usual realization of the quark mass term. The same overall picture 
was fomid. However, because the method relied on the symmetry properties of the Lagrangian, only the properties 
of the pseudoscalar states and the strange scalar states could be studied. In the present paper, we shall initiate a 
much more systematic investigation. We first study precisely how the general results get constrained when a specific 
choice of invariant interaction terms is made. We introduce a scheme for ordering all the non-derivative terms of the 
Lagrangian according to their likely importance. This enables us to select a limited number of leading order terms 
in a meaningful way as well as to provide the framework for possible higher order extensions. At leading order, and 
with the extra simplification of zero quark masses, all our results were determined analytically, without any need for 
a numerical fitting procedure. There are essentially only four main input parameters and only one of them has a non 
negligible experimental error. We do our calculations for all allowable values of this parameter (m[7r(1300)]) and also 
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take into account the small experimental error on another of the three parameters. The results obtained dramatically 
predict the existence of a very low mass scalar isosinglet state. Especially, the puzzle concerning the coexistence of 
lighter (mainly) two quark pseudoscalars with lighter (largely) four quark scalars is clearly seen to be solved. 

A brief review of the model and the relevant notation is presented in Sec. II. Section III shows the great simpli- 
fications obtained by going to massless QCD and also gives our notations appropriate to the flavor SU(3) invariant 
sitTiation in this limit. General results, valid for any choice of terms in the invariant potential, are also presented in 
this section. Section IV gives a systematic procedure for deciding which terms are most important in the model. It 
mainly contains the worked out model using the leading terms in this scheme. A numerical analysis is presented and 
the masses of the two SU(3) singlet scalar states of the model are predicted. The two and four quark contents for each 
state of the model are displayed. Sections V, VI and VII are devoted to proving, for any choice of invariant potential, 
the current algebra theorem for the scattering of massless pions. Discussion and conclusions are given in Sec. VIII. 



II. BRIEF REVIEW OF MODEL 



The fields of our "toy" model consist of a 3 x 3 matrix chiral nonet field, M which represents qq type states as well 
as a 3 X 3 matrix chiral nonet field, M' which represents qqqq type states. They have the decompositions into scalar 
and pseudoscalar pieces: 

M = S + i(t), 

M' = S' + icf)'. (1) 
They behave under "left handed" and "right handed" unitary unimodular (ie SU(3)lX SU(3)r) transformations as 

M' UlM'uI^. (2) 

However, under the U(1)a transformation which acts at the quark level as qaL — * e^^^aL-, QaR ^ ^~^'^<laRj the two 
fields behave differently: 

M' e-^'-^M'. (3) 

We will be interested in the situation where non-zero vacuum values of the diagonal components of S and S' may 
exist. These will be denoted by, 

{Si) = aJl {S'^) = f3Jl (4) 

In the iso-spin invariant limit, ax = a2 and /3i = P2 while in the SU(3) invariant limit, ai = a2 = as and (3i = (^2 = Ps- 
The general Lagrangian density which defines our model is 

£ = -^TV (a^Ma^Mt) - ItV (a^M'S^M't) - Fo (M, M') - VsB, (5) 

where Vo(M, M') stands for a general function made from SU(3)lX SU(3)r (but not necessarily U(1)a) invariants 
formed out of M and M' . The last term, Vsb, stands for chiral symmetry breaking terms which transform in the same 
way as the quark mass terms in the fundamental QCD Lagrangian. In the present paper we shall, in later sections, 
specialize to the zero quark mass limit by taking Vsb = 0. Not only does this make the formalism much simpler but 
it is well known that, due to the spontaneous breakdown of chiral symmetry, the main mechanism of physical hadron 
mass generation is already accounted for. This is convenient for disentangling the general properties of each multiplet 
from the imcertainty as to which of the many possible mass type terms in the cfffcctive Lagrangian to include. We 
record the behaviors of the fields under infinitesimal transformations. Let us write the infinitesimal vector (L+R) and 
axial vector (L-R) transformations of (j) and S as, 

5v(t> = [EvA]-, 5a4> = -i[EA-,S\+, 

SvS = [Ev,S], 5 AS = i[EA, <!>]+. (6) 
Here, unitarity demands that the infinitesimal matrices obey, 

El = -Ev, e\ = -Ea. (7) 
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If we demand that the transformations be unimodular, so that the U(1)a transformation is not included (the U(l)v 
transformation is trivial for mesons), we should also impose Tr(i?yi) — 0. However we will not do this so the effects 
of U(1)a will also be included. The transformation properties of the qqqq type fields are: 



Sy(j)' = [Ev,(j)']- 

6vS' = [Ev,S'] 



SacIj' = ~i[EA, S']+ + 2iS'Tt{Ea), 



SaS' = i[EA, ( 



2i0'Tr(£;A). 



(8) 



The extra terms for the axial transformations reflect the different U(1)a transformation properties of M and M' . 

We will employ two complementary approaches to make predictions. One approach will be to study generating 
equations for tree level vertices. These are like Ward identities and follow for any choice of Vq{M,M') in Eq. ([S]). 
These predictions are consistent with but will not give all possible predictions which would arise if one considered, as 
a second approach, making a specific choice of terms in Vb(M, M'). 

The method of treatment, as used earlier [i^ to discuss the model containing only the field M, is based on two 
generating equations which reflect the invariance of Vp under vector and axial vector transformations. Differentiating 
them once, relates two point vertices (masses) with one point vertices. Differentiating them twice relates three point 
vertices (trilinear couplings) with masses and so on. Under the infinitesimal vector and axial vector transformations 
we have, 



^v^^ = Tr ( ^5v<i^ + -Q§5vS ) + (0, S) ^ (</>', S') = 0, 
' BV BV 

SaV^ = Tr ( -^5A<ly + -q^SaS ) + (0, S) ^ (0', S') - -C^, 



(9) 



wherein the non-zero value of the axial variation equation reflects the presence in Vq of any terms which are not 
invariant under U(1)a; these terms will provide mass to the ?7'(958) meson. In (46j . terms of this type were represented 
by any function of the chiral SU(3), but not U(1)a, invariant det(M) plus its hermitian conjugate. After QCD, 't 
Hooft found [131 that such a form would arise from instanton effects. If one wishes to model the U(1)a anomaly 
equation of QCD in the single M model the suggested form is: 



-C3 



In 



/ detA/ \ 
V detMt J 



(10) 



where C3 is a numerical parameter. In the present M — M' model this form is not unique and the most plausible 
modification [4§| is to replace ln(^^^^) by 



71 In 



/ det(M) \ 
Uet(Mt) J 



/Tr(MM't) 
(^-^^)^HTr(M'Mt) 



(11) 



where 71 is a dimensionless parameter. Using Eqs. ^ and ([5]) as well as the arbitrariness of the variations Ey and 
Ea yields the matrix generating equations, 



' d(i) 

' 8S 



' dVo 
8S 



+ (0,5)^(0', 5') = 0, 
+ (0,5)->(0',5') = l 



2Tr 



,dVo ^,dVo\ ^ ., /dctM\ 



(12) 



where the form of Eq. pH)) was used. In addition, the replacement, Eg. pTj) should be borne in mind. To get 
constraints on the particle masses we will differentiate these equations once with respect to each of the four matrix 
fields: 4>, (f>\ 5, S' and evaluate the equations in the ground state. Thus we also need the "minimum" condition. 



/dVo 

\ 8S 



SB 



8S 



= 0, 



m 

8S' 



/ dVsB 
\ dS' 



= 0. 



(13) 



In ref [4J| we considered the canonical term, Vsb = — 2Tr(yl5) as an effective representation of the fundamental quark 
mass terms; A is a diagonal matrix with entries proportional to the three quark masses. Next, let us differentiate 
successively the axial vector generating equation with respect to </> and to (j)' . It is neater to write the results first for 
the case when fields with different upper and lower tensor indices are involved: 
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{Ua + at) 



{Pa + 1%) 



+ [Pa + f3b) 



= 



(14) 



Next, let us write the corresponding equations for the case when the upper and lower tensor indices on each field 
are the same. 



at 



-Pb 



a 

E 

g 



d<i>id^'l 



8C3 
Ota ' 



(15) 



Note that the axial generating equation provides information on the masses of all the pseudoscalars. Further differ- 
entiations will relate a large number of trilinear and quadrilinear coupling constants to the meson masses and to the 
quark mass coefficients, Aa- 

To fully characterize the system we will also require some knowledge of the axial vector and vector currents (46j 
obtained by Noether's method: 



where the dots stand for terms bilinear in the fields. 

It will be helpful to briefly review the treatment of the it-it' system as given in section IV of ref. [ii 
the abbreviations 

ai 



(16) 
Introduce 



ai ' 



(17) 



Here we have introduced the total potential V = Vq + Vsb- Substituting a = 1, 6 = 2 into both of Eqs. (fT4)) enables 
us to write the (non-diagonal) matrix of squared tt and tt' masses as: 



(18) 



We see that a;^ would be the squared pion mass in the single AI model and j/tt represents the squared mass of the 
"bare" tt' . The transformation between the diagonal fields (say 7r+ and 7r'+) and the original pion fields is defined as: 



r'+ 



cos ( 
sin^ 



— sm( 
cos 0. 



(19) 



which also defines the transformation matrix, R. The explicit diagonalization gives an expression for the mixing angle 



tan(26'^) = 



y7r(l - - X^' 



(20) 



The mixing angle, can be connected to the experimentally known value of the pion decay constant. Substituting 
the expressions from Eq. for (f)\ and 0'^ in terms of the physical fields 7r+ and 7r'+ into Eq. yields. 



Ft, = (ai + 02) cos6'^ - (/3i -I- /32) sin^^, 
F-n- = (ai -I- 02) sin6'^ -I- (/3i -I- /32) cos6'^. 



(21) 
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III. SIMPLIFICATION FOR ZERO QUARK MASSES 

The zero quark mass limit is gotten by taking Vsb — 0. We assume that the original SU(3)lX SU(3)r symmetry 
is spontaneously broken to SU(3)v rather than some smaller subgroup. The vacuum expectation values of the scalar 
fields simplify to: 

ai = "2 = as = a, Pi = 02 = Ps = P- (22) 

The mass spectrum also simplifies a lot. When quark masses are included in the isotopic spin invariant approximation 
there are 16 different masses. However in the zero quark mass limit there are only 8 different masses. These describe 
the four systems of degenerate SU(3) octet or SU(3) singlet fields: 

(0,0'); (00, 0o), 

(^,^0, (^o,^o)- (23) 

Here the hat stands for the eight members of the appropriate octets. The fields of each system can mix with each 
other but not with the fields of any other system. In addition to 8 different masses there will be four different mixing 
angles describing four orthogonal 2x2 matrices. The conventions are the same as in Eq. (fT9|) so that 9^ now describes 
the mixings of the two pseudoscalar octets. Note that if only isotopic spin invariance were present, the isotopic spin 
zero fields of each parity would be characterized by a 4x4 mixing matrix with 6 angle parameters (See Eq.(64) of |44| 
for example). Notice also that the tt — tt', K — K' and ?78 — ?7g mixings, for example, are all described by the same 
mixing parameter 6^^. 

We next discuss the notations for resolving the nonets into SU(3) octets and singlets. Matrix notation is sometimes 
convenient; we use the convention cj)'^ — > (j)ab- The properly normalized singlet states are: 

00 = ^Tr(0), q^'o = ^Tr(0'), 

^0 = ^Tr(5), S'o = ^MS')- (24) 



Then we have the matrix decompositions: 



S = S+^Sol, S' = S' + ^S'l, (25) 

VS \/3 



wherein (p, (j)', S and S' are all 3x3 traceless matrices. The singlet scalar fields may be further decomposed as: 

So^VSa + So, S'o = VS(3 + S'q. (26) 

Here Sq and S'q are the fiuctuation fields around the true ground state of the model. 

Setting = 0, corresponding to zero quark masses, simplifies Eq. (PO)) for the tt-tt' mixing angle to: 

, o/j -2z7r _ 2tan6l^ 



This immediately yields: 



tan6l^ = --. (28) 
a 



Substituting this into Eq. (PTjl yields the simple results: 



= 2Va2+/32, 

F^, = 0. (29) 

One may note, for comparison, from Table 4 in [45| that F^^i and also Fx' are not exactly zero in the presence of 
non zero quark masses, although they are very heavily suppressed. This feature suggests the essential reliability of 
the zero quark mass limit. 
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Next consider the pseudoscalar octets, 
just the two / = /3 = 1 fields, 0f and 



(f> and (/)' in the model. Because of SU(3) symmetry it is sufficient to give 
Their mixing matrix, Eq. (jlSp becomes in the limit of zero quark masses: 



(Mi) 



-z^ 1 



-13/a 



-(3/a 
1 



(30) 



It is easy to see that this matrix has zero determinant and to identify the usual (but zero mass) pseudoscalar pion as 

2 



(31) 



where =131 MeV. The physical massive pion "excitation" is clearly 7r'+ = ■^(— /90^ + acj)' i) and has a squared 
mass, m?'{iT') = yTr(l + 0^ /cP'). We notice that, just from our general treatment, the tt — it' system can be described 
by the three parameters a, (3 and j/tt- However, there are only two physical quantities, and rnPiT:'), to compare 
with. Thus the mixing angle between the usual and the "excited" pseudoscalar octet states is not predicted in 
general. In order to predict this interesting quantity we have to specify our choice of chiral invariant terms in the 
potential, V . A similar situation will be seen to hold for trilinear and quadrilinear coupling constants involving the 
physical pseudoscalars. There are many constraints just from chiral symmetry but a complete (though clearly model 
dependent) description will depend on the particular choice of terms in the potential. 

It is also amusing to look at the 00 — 0'o sector in the zero quark mass limit. There is a rather drastic simplification 
since the introduction of quark masses results in additional mixing with the isoscalar members of the corresponding 
octets. That requires a six parameter 4x4 transformation matrix rather than the single parameter 2x2 matrix we 
now will get. Using the formula, 

\ I *^ *^ \ 

(32) 



d ^ 1 


r d 


d 


d 


dcbo V3 




d4>l ^ 


d<t>l 



in both of Eqs. (|15p . we end up with the pre-diagonal 0o ~ 0'o mass squared matrix: 



(Mo! 



- ^^^^^ 



8c3(l-7i)(27i + l) 
3^2 -202/0 + ^ 

^ „ , 8c3(l-7l)(27l + l) „ 8c3(l-7l)'' 
-ZoVa H ^ 2/0 3^2 



Here zq = — 2/3/a and 



yo 



_&^V_ 



The mixing angle, O^j, is defined by the convention: 



0Op 




00 




>0P. 


.0'o. 





cos 00 — sin 00 
sin 9a cos 9q 



(33) 



(34) 



(35) 



In the limit where C3, defined in Eq. (jlOp. vanishes it is seen that the determinant of the mass squared matrix in Eq. 
(j33p vanishes. This is understandable since C3 multiplies the terms which violate U(1)a symmetry and a zero mass 
singlet pseudoscalar boson must exist since the symmetry is broken spontaneously. In this limit the mixing angle is 
related to the pseudoscalar octet one by tan^o = —2 tan 6*^. It should be remarked that the effect of non zero C3 is 
actually quite large so the limit where it vanishes is mainly of academic interest. 

Of course, the U(1)a transformation is relevant in setting up this model since, as seen from Eq. ([3]), it distinguishes 
the two quark fields from the four quark fields. We shall consider here, models in which the terms multiplied by C3 
are the only ones which violate U(1)a symmetry. In that case the divergence of the axial current in the model exactly 
mocks up the QCD axial anomaly at tree level. Alternatively, a term like det(M) + h.c. could be used with similar 
results; such a term does not however mock up the U(1)a anomaly equation. 

We have seen that quite a lot of information about the pseudoscalar particle masses and mixings follows just from 
the axial generating equations, reflecting the spontaneous breakdown of the octet axial symmetries. On the other 
hand, Eq.(31) of ref.[44[ shows that, in the case where spontaneous breakdown preserves the SU(3) invariance of the 
vacuum, there will be no such model independent information about the masses and mixings of the scalars. To find 
that information, one must make models with specific choices of the invariant terms. In preparation we give notations 
for the scalar mass and transformation matrices, analogous to those we adopted for the pseudoscalars, in the case 
where quark masses are absent and the vacuum is assumed to be SU(3)v invariant. The pre-diagonal 2x2 matrix 
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for the 7 = 73 = 1 scalar meson squared masses is denoted (X^) and the mass diagonal fields, a"*" and a'+ are related 
to the non-diagonal ones by: 



1'+ 



= L- 



S'l 



a 

— sin ipa 

Sinipa cos Ipa 



si 

S'l 



(36) 



This is sufficient to describe the mixing of all the scalar octet particles with corresponding SU(3) quantum numbers. 
For the So — S'q mixing, we define the prediagonal squared mass matrix to be {Xq ) while the mass diagonal fields 
and are defined by: 



'Op 



Op 



1 


5o 




cos V'o 


— sin V'o 




'So' 








sin V'o 


cos ipo 







(37) 



IV. MODEL FOR MASSES AND MIXINGS 



As just discussed, it is necessary to make a specific choice of terms in the SU(3)lX SU(3)i^ invariant potential Vq 
in order to be able to predict all physical properties of the system. This is a non trivial issue since, for example, if we 
restrict Vq to be renormalizable, there are twenty one terms (4^ with this symmetry. We will adopt two criteria for 
which terms to include. First we list the six SU(3)lX SU(3)i^ invariant terms which satisfy these criteria and seem 
the most reasonable for an initial treatment: 

Vo^ - C2 Tr(MMt) + cl Ty{MM'^MM'^) 

+ d2 TiiM'M'^) + {eabc^''^^ M^M^M'f^ + h.c) 



C3 



, /detM\ . /Tr(MM't) 
^^^Hd^J^^ ''^^ l Tr(M-Mt) 



(38) 



All the terms except the last two have been chosen to also possess the U(1)a invariance. Those terms are clearly 
non-renormalizable and violate U(1)a invariance in a special way. They have, as previously discussed [see Eq. ljlip ]. the 
correct U(1)a property so that the resulting Lagrangian can exactly mock up the U(1)a anomaly of QCD. Of course, 
we are using the effective Lagrangian at tree level and renormalizability is not an issue at this level. Renormalizable 
terms of the instanton determinant type and the type Tr(MM't) + h.c. could be used instead with not much change 
in the result. However, the role that the U(1)a transformation is playing in distinguishing "four quark" from "two 
quark" effective fields suggests that we try to reproduce as much as possible of the behavior of QCD under axial 
U(1)a. The In terms chosen also have the convenient feature that they confine the U(1)a violating effects to the 
SU(3) singlet pseudoscalar sector of the model. The first four terms were chosen from the twelve renormalizable 
and U(1)a invariant ones in the formula, Eq. (Al) of [31 (please see also Appendix A of the present paper) by 
imposing the criterion that effective vertices describing the smallest numbers of quarks plus antiquarks be retained. 
This quantity, representing the total number of fermion lines at each effective vertex can be written as, 

7V = 2?i + 4n', (39) 

where n is the number of times M or appears in each term while n' is the number of times M' or M'^ appears 
in each term.. Thus, the C2 term has N=4 while the cf, c?2 and eg terms each have N=8. For simplicity, we have 
neglected the 7V=8 term, C4[Tr(MMt)]^ which is suppressed, in the single M model, by the quark line rule. It may be 
noted that the quantities det(M) and Tr(MM't) which enter into those two terms which saturate the U(1)a anomaly 
have A'^=6. On the other hand, the terms in 

64 Tr(MA'7tM'A/'t) + el Tt{MM'^ M' M'') (40) 

each represent twelve quarks plus antiquarks at the same vertex and will not be included at the present stage. Similarly, 
the term c?f Tr(M'M'^M'M'^) representing sixteen quarks and antiquarks will not be included. In the future, U(1)a 
invariant terms with higher values of N may be used to systematically improve the approximation as well as U(1)a 
violating operators with higher values of N which may be inserted into an obvious generalization of Eq. pip . The 
minimum equations for this potential are: 

-2a (-C2 + 2c2a2+4e^/3) -0, (41) 
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dVo 



2{d2P + 2e^a^) = 0. 



(42) 



Notice that a is an overall factor in Eq. (|4T|) so that, in addition to the physical spontaneous breakdown solution 
where a ^ there is a solution with a = 0. On the other hand, (3 is not an overall factor of Eq. (I42|l and it is easy 
to see that /3 is necessarily non-zero in the physical situation where a is non-zero. The minimum equations clearly 
eliminate two parameters from the model. 

Next, we shall give the matrix elements of the four squared mass mixing matrices based on the use of the specific 
potential of Eq. (|38p . First consider the matrix describing any of the eight degenerate 0~ quark-antiquark fields mixing 
with their corresponding four quark partners. Without using the minimum equations, one obtains: 



2 (-C2 + 2 ct + 2 eg /3) 4 eg a 

2d2 



4 eg a 



This corresponds to the general form given in Eg. pop when we identify, y^r = 2d2 and z^^ 
= (3 /a. 

The matrix describing the mixing of the two pseudoscalar singlets is similarly written as: 



(43) 

-2ae1/d2. Note that 



(Mo) 



-2 (c2-2c4-a2+4eg/?)-«23(|^ 



2d2 



8c3(l-7l)(27l+l) 

8c3(l-7l)'' 
3p 



(44) 



This corrsponds to the general form given in Eq. (j30p when we identify, j/o ~ 2^2 and zq — —2f3/a — 4ega/d2- 
For the mixing matrix of the octet scalars, the specific potential of Eq. ((38|) directly gives: 



2 (-C2 + 6cta2-2eg/3) 



-4ae3 



-4Q;e3 
2d2 



Finally the squared mass mixing matrix for the singlet scalars is similarly obtained as: 

2 (-C2 



6 c| a2 + 4 ef /3) 8aef 

2d2 



8aeg 



(45) 



(46) 



Now let us consider the comparison of this model with experiment. To start with there are 8 parameters (a, /3, C2, 
•^2, C4, eg, C3 and 71). These can be reduced to six by use of the two minimum equations just given. We note that 
the parameters C3 and 71), associated with modeling the U(1)a anomaly, do not contribute to either the minimum 
equations or to the mass matrices of the particles which are not 0~ singlets. Thus it is convenient to first determine 
the other four independent parameters. As the corresponding four experimental inputs [l| we take the non-strange 
quantities: 



TO(0'''octet) 
m(0+octet') 
m(0~octet') 



m[ao(980)] = 984.7 ± 1.2 MeV 
m[ao(1450)] = 1474 ± 19 MeV 
m[7r(1300)] = 1300 ± 100 MeV 
131 MeV 



(47) 



Evidently, a large experimental uncertainty appears in the mass of 7r(1300); we shall initially take the other masses 
as fixed at their central values and vary this mass in the indicated range. As shown in Eq. (jBip in Appendix B, it 
is straightforward to determine the four independent parameters in terms of these masses. There is a complication 
which must be taken into account; from studying the predicted masses of the O''" SU(3) singlet states one finds that 
the positivity of the eigenvalues of their squared mass matrix, Eq. (|46p is only satisfied when, 

(48) 



[7r(1300)] < 1302 MeV. 



Further restrictions on the allowed range of r7i[7r(1300)] will arise when we calculate the masses of the SU(3) singlet 
states. Before that we mention the two predicted masses for the O'^ SU(3) singlet states; as m[7r(1300)] varies from 
1200 to 1300 MeV, 



m(0''' singlet) 
m(0"'" singlet') 



510 -> 28 (410) MeV, 
1506 ^ 1555 (1520) MeV. 



(49) 
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m [ltd 300)] (GeV) 

FIG. 1: The predictions for the masses of the two SU(3) singlet scalars vs. m[7r(1300)]. The error bars give the effect of the 
uncertainty in the ao(1450) mass. 

The predictions in parentheses correspond to the hkely additional constraints from the positivity of the 0^ SU(3) 
singlet states. Plots are shown in Fig. [T] 

Clearly, the most dramatic feature is the very low mass of the lighter SU(3) singlet scalar meson. Of course, one 
expects the addition of quark mass type terms to modify the details somewhat. On the other hand, there are a number 
of allowed different quark mass terms so it is notable that the characteristic very light mass scalar exists apart from 
the ambiguity in choice of the quark mass terms. 

The four independent parameters which appear in the Lagrangian (c2, d2, C4, 63) are shown, as functions of 
TO[7r(1300)], in Fig. O The vacuum expectation values a and /3 of the two and four quark scalar fields are similarly 
shown in Fig. [31 It is seen that /3 and a are each insensitive to varying m[7r(1300)] and their ratio is about 0.40. 

To calculate the masses of the SU(3) singlet pseudoscalars we must diagonalize Eq. (|55)) with the specific choices of 
parameters yo = 2d2 and zq = 4:e^a/d2 corresponding to the potential of Ea. (|38p . This enables us to fit in principle, 
for any choice of m[7r(1300)], the two parameters C3 and 71 in terms of the experimental masses of 77(958) and one 
of the candidates 77(1295), 7^(1405), 7/(1475) and 7;(1760). The specific formulas are given as Eqs. (|B2p and (jB3[) in 
Appendix B. However, as mentioned above, the positivity of the eigenvalues of the matrix (Mq) imposes additional 
constraints on the choice of 77T,[7r(1300)] in Eq. ipS)) . This appears in solving for 71 using the quadratic equation (|B2|) 
and requiring its discriminant to be positive. In Fig. [H the discriminants are shown as functions of 77i[7r(1300)] for 
each of the four possible candidates for the heavier 0~ SU(3) singlet. This clearly shows that the two lowest mass 
candidates have negative discriminants and can be ruled out according to our criterion. The perhaps most likely 
candidate 7/(1475) [this case will be denoted scenario 1] has a positive discriminant for 77i[7r(1300)] less than about 
1.23 GeV. This leads to the modified allowed ranges for the 0+ singlet states, shown in parentheses in Eg. (1^5]) . There 
is no restriction on the heaviest candidate, 77(1760) [this case will be denoted scenario 2]. 

Since Eq. (jB2p is a quadratic equation for 71, one expects that there may be two physical solutions for 71. This 
turns out to be the case. In FigOwe show plots of 71 as a function of 777,[7r(1300)] for each of the scenarios mentioned 
above. The quantity C3 is given in Ea. (jB3p and is seen in Fig|S]to be single valued in its dependence on 777,[7r(1300)]. 

It is very interesting to see what the model has to say about the four quark percentages of the particles it describes. 
The percentages for the pion, the lighter 0+ singlet and the ao(980) are displayed in Figl7]as functions of the precise 
value of the input parameter 77z[7r(1300)]. The pion four quark content (equal to 100 sin^0^) is seen to be about 17 
percent. Of course the heavier pion would have about an 83 percent four quark content. On the other hand, the octet 
scalar states present a reversed picture: the ao(980) has a large four quark content while the ao(1450) has a smaller 
four quark content. The very light and the rather heavy 0''" singlets are about maximally mixed, having roughly equal 
contributions from the 4 quark and 2 quark components. 

In FigjS] the four quark percentages of the 0~ SU(3) singlets are shown for both scenarios. The perhaps more 
plausible scenario takes 7/(1475) as the heavy 0~ singlet state. In this case we see that for the solution with smaller 
71, the four quark content of the familiar 7/(958) is about 25 percent while for the solution with larger 71, the four 
quark content of 77(958) is about 55 percent. Thus the smaller 71 solution seems more plausible physically. In the 
case where the 77(1760) is identified as the heavier partner of the 77(958) the smaller 71 solution yields an 7/(958) with 
a four quark content of about 7 percent while the larger 71 solution yields an 7/(958) with a four quark content of 



1.2 1.22 1.24 1.26 1.2S 1.3 1.2 1.22 1.24 1.26 1.28 1.3 

m [1(1300)1 (GeV) m [lt(1300)] (GeV) 




m [It (1300)] (GeV) ra [1(1300)1 (GeV) 



FIG. 2: Starting from the upper left and proceeding cloclcwise: C2 vs. m[7r(1300)], d2 vs. m[7r(1300)], ef vs. m[7r(1300)] and cf 
vs. m[7r(1300)]. The range of m[7r(1300)] corresponds to the restrictions imposed by the positivity of the scalar SU(3) singlet 
masses. The error bars give the effect of the uncertainty in the ao(1450) mass. 



about 82 percent. 

Values of all the model parameters as well as numerical values of the mixing matrices, for a typical choice of 
TO[7r(1300)], are hsted at the end of Appendix B. 



V. THREE POINT VERTICES 



The three point vertices are useful for calculating the widths of the various mesons and also for the calculation of 
meson-meson scattering. These can be calculated for a specific model, like the one with the choice of terms given in the 
previous section, by straightforward differentiation. However, one may also obtain model independent (in the sense 
of being independent of the choice of invariant terms in Vq) information about these from the generating equation. 
We shall do that here, specializing to the scalar-pseudoscalar-pseudoscalar vertices needed for pion pion scattering. 
These are obtained by succesively differentiating the two equations in Eq. (jl2p with respect to one scalar field and one 
pseudoscalar field. First we introduce the notations: 

d^Vo \ I a^^o 



^1 \ PiA.lp,A.'2-F\Q. I ^1 



02'J(p\dSii I \ d(j)\dcf)\dS% 

a^Fo \ / a^vb 



''^ ^ ''h\d(j)ldsu \d^\d^ldS'i 
^ _ I d^V^ \ _ / a^vb 

"""^ \d(l)ld(l)[^dSo / '^^ \d(f>^2d(f>'ldSi 



m |II(1300)J (GeV) 



m [1(1300)1 (GeV) 



FIG. 3: Dependences of the two quark vacuum value a (left) and the four quark vacuum value (5 (right) on the choice of 
m[7r(1300)]. The error bars give the effect of the uncertainty in the ao(1450) mass. 




FIG. 4: The discriminant for Ea. (|B2|l vs. m[7r(1300)]. The curves from bottom to top respectively represent the choices for 
the heavier 0~ SU(3) singlet to be ry(1295), r;(1405), j7(1475) and 7^(1760). The error bars give the effect of the uncertainty in 
the ao(1450) mass. 



d^Vo \ I a^vb 



Note that 5*0 was defined in Eg. ((24)) while Sg, for example, is the isoscalar member of the SU(3) octet defined as: 

S^ = ^^{S\-rSl~1Sl). (51) 

Now using the generating equations as just discussed, we obtain the following relations connecting the trilinear 
coupling constants r,; with corresponding mass squared matrices for the Sq-Sq and the tt-tt' systems. 



ari + jSrs 



/d^Vo\ 1 / d'Vo 
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1.22 

111 [JI(1300)] (GeV) 



FIG. 5: 7i vs. m[7r(1300)]. The top and bottom curves correspond to choosing the »7(1760) as the heavier 0~ SU(3) singlet 
while the middle two curves correspond to choosing the 77(1475) as the heavier 0~ SU(3) singlet. Note that for each scenario, 
the two curves are associated with different solutions of the quadratic equation IB2I for 71. The error bars give the effect of the 
uncertainty in the ao(1450) mass. 




m [71(1300)1 (GeV) 



FIG. 6: C3 in units of GeV^ vs. m[7r(1300)]. The upper curve corresponds to the scenario where the heavier SU(3) singlet 
is identified with the 77(1475) while the lower curve corresponds to 77(1760) as the heavier 0~ SU(3) singlet. 



ar2 



(3r4 
Pre 



1 



1 / d'-Vo 
V3\dSodS'o 
1 / d^Vo 



d'Vo 



VS \ dct>ld<j>'l 
J_/_d^Vo 
V3 



V3\dS'l/ V^\d(P'ld(l)'l 



(52) 



Similar equations are obtained for the qi trilinear couplings and the mass squared matrices for the Sg-S'g systems: 



aqi 



aq2 + Pq4 = 



1 /d^Vo 
V6 \ dSi 
1 / d^Vo 

V6 \ dSsdS's I \ 50^(90'i 



1 / a^^o 

1 IJ^ 



otq-i 



_ 1 / g^Fo \ 1 
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1.24 

m [71(1 300)] (GeV) 



FIG. 7; Plot of the four quark percentages of various particles in the model as functions of the undetermined input parameter, 
m[7r(1300)]. Starting from the bottom and going up, the curves respectively show the four quark percentages of the pion, the 
0^ singlet, and the ao(980). The error bars give the effect of the uncertainty in the ao(1450) mass. 




Ill [11(1300)] (GeV) 



FIG. 8: Plot of the four quark percentages of the »7(958) as functions of the undetermined input parameter, m[7r(1300)] for two 
scenarios. The top and bottom curves correspond to choosing the 77(1760) as the heavier 0~ SU(3) singlet while the middle 
two curves correspond to choosing the 77(1475) as the heavier 0~ SU(3) singlet. Note that for each scenario, the two curves are 
associated with different solutions of the quadratic equation IB2I for 71 . The error bars give the effect of the uncertainty in the 
ao(1450) mass. 



aq4 + Pqe 



V6 \ dS'\ 



V6 



(53) 



Eq. ((52l) and Eg. ([53]) relate eight different linear combinations of the three point vertices to two point vertices for 
the fields of pure qq and pure qqqq types. Since there are twelve a priori unknown three point vertices according to 
Eg. ([50)1 ■ it is clear that there is, in general, not enough information available to determine all the three point vertices 
in terms of the two point ones. However, we will see that the available relations are sufficient to prove the desired 
low energy theorem. To relate the quantities in Egs. ([52)) and ([53)1 to quantities pertaining to mass eigenstates we 
introduce an index notation to distinguish unprimed from primed fields; for example: 



{<P\)2- 



(54) 



With this notation, which we apply to all fields of the model, the coupling constant of the Goldstone boson pions to 
the mass diagonal SU(3) singlet scalars may be compactly written as: 



ffOD 



\ dTr+dTr-d{Sop)D 



— ^ (-Rir)Al(-Rir)si(io) 



CD 



A,B,C 



d{(p\)Ad{^\)Bd{So)c 



(55) 
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The transformation matrix elements, {Rtt)ab and {Lo)ab may be read from Eq. (fT9)l and Eg. ([37)1 . Note that the 
capital Latin subscripts take on the values 1 and 2 as in Ea. (|54p above. There is a similar equation involving the 
S'g-S'g scalars which yields the physical coupling constant of two Goldstone pions with 5*8, gso- 



gsD 



\ / QSy \ 

dir+dTT d{Ssp)D/ \d{(l)i)Ad{cl)^)Bd{Ss)c / 



Here, L is the transformation matrix defined in Eg. pop . Using the compact form of Eq.(|55p. one may compactly 
express the comparison of Eg. ([5^ with Eg. ([SO)) as: 

y2i^^^)iB ( c,(<2\ fUJA\ afo \ \ = i^o)AH - (M^)Aif. (57) 

2 ^ \d{(t)f)Ad{(l)^)Bd[So)H / 

(M^) is given in Eg. ([50)1 and (Xq) is the model independent version of Eq. (P5|) . Similarly, 

^ " y2i^^^)iB ( p,,,2\ f^(Q \ ') = i^a)AH - {MI)aH- (58) 

2 ^ \d{(t)i)Ad{(t)'2)Bd{Ss)H / 

Here [X"^) is the model independent version of Eg. ([45]) . Note that according to our conventions the nondiagonal and 
diagonal (hatted) squared mass matrices are related as: 



Y,{R-^)ab{MI)bc{R.)cd = {MDad, Y.^Ro^^ab{M^)bc{Ro)cd = {M^)ad, 

B,C B,C 

^{La^)AB{Xj)Bc{La)cD = {XI)aD, ^{L^^)ab{Xq)bc{Lo)cD = {Xg)AD, (59) 
B,C B,C 

VI. LOW ENERGY PION SCATTERING 

There are two reasons for next discussing the pi-pi scattering in this model. First, since the iso-singlet scalar 
resonances above are being considered at tree level, one expects, as can be seen in the single M model also discussed 
in [3l[ and at the two flavor level in Q , that unitarity corrections for the scattering amplitudes will alter their masses 
and widths. Second, since the pion looks unconventional in this model (having a non-neglegible four quark component) 
one might worry that the fairly precise "current algebra" formula for the near to threshold scattering amplitude might 
acquire unacceptably large corrections. 

Of course, for computing the near threshold pion pion scattering, it is well known that the use of a nonlinear sigma 
model is more convenient. However, we are also interested in unitarizing the model in the resonance region where the 
nonlinear model, which can be obtained by integrating out the resonances, is clearly not applicable. 

The invariant pion pion scattering amplitude for 'Ki{pi)+T: j{p2)^ '^kiP3)+'^i{p4:) is decomposed as: 

6ijSkiA{s, t, u) + 5ik6]iA{t, s, u) + 6ii5.jkA{u, t, s), (60) 

where s, t and u are the usual Mandelstam variables. Note that the phase of the above amplitude simply corresponds 
to taking the matrix element of the Lagrangian density for a four pion contact interaction. The / = 0, / = 1 and 
1 = 2 amplitudes correspond to the projections: 

T°{s,t,u) = 3A{s,t,u) + A{t,s,u) + A{u,t,s), 
T^{s,t,u) = A{t, s,u) — A{u,t, s), 

T^{s,t,u) = A{t,s,u) + A{u,t,s). (61) 

It is straightforward to calculate A(s, t, u) using the three point vertices for two massless pions coupling to a physical 
scalar (See Eqs. ([55]) and ([56)) ) as well as the four point coupling constant, g for four massless pions: 

> (62) 



dn+dn dTT+dn 



The result is simply: 

JsD ^ 9o_ 

[XDdb - s {XI)dd - s 



A{s, t,u) = -'- + J2( j^^^ + j^^^ ] ■ (63) 
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Note that the sum goes over the two SU(3) singlet scalars as well as the two iso-singlet scalars belonging to SU(3) 
octets. We are presently interested in the threshold region (near s=0 for massless pions) so we expand this formula 
to first order in s: 

A(s,t,u)^~^+(^^ + ^^]+.s(^^ + V (64) 

2 \iX^^)oD iXi)nDj \m)DD? [iXi)DD?J 

In this equation the summation over D has not been explicitly written and the summation over repeated indices is 
to be assumed; note that the quantity (X^)dd, for example, is a single number indexed by D. Observe that the four 
point vertex does not contribute to the terms linear in s. Let us then evaluate the s term first. Begin by substituting 
Ea. (|57p into Eq. (|55p and noticing that the term {M^)ah makes zero contribution since that piece can be manipulated, 
using Eg. ([55]) . to be proportional to the zero masses of the physical Goldstone bosons. The physical trilinear coupling 
constant is next obtained as goo = {R-r)ai{Xq)ah{Lq)hd- Then the quantity appearing in Eg. ([64]) can be 

evaluated as: 



9aD 



HD 



{Rtt)ci{Xq)ck{Lo) 



KD 



)DD\ 



3F2 
4 

m 



;i^o ^)dk{xI)kc{Lo)cf{Lq ^)Fj{R-n).n 



DD\ 



{RTr^)iG{La)GE{Xl) 



ED 



DD 



-{xI)df{Lq^)fj{R7t)ji - 



Similarly, 



91d 



4 



The s dependent part of the scattering amplitude near threshold finally takes the simple form: 

A{s,t,u) = — . 



(65) 



(66) 



(67) 



This may be recognized as the usual current algebra formula [50[ in the case where the pion mass is set to zero. 
We will complete its derivation in the next section, where it will be shown that the s independent terms in Eg. (j64p 
cancel each other. It should be remarked that the present derivation holds for any choice of chiral invariant terms in 
Vo, not necessarily just for the leading terms in Eg. ([55)1 . 

Of course, the current algebra result is just the first term in an expansion in powers of s. The higher terms will 
have the structure of a geometric series: 



A(s, t, u) 



9 
m 



(68) 



wherein we have amalgmated all four scalars as the nii and their corresponding coupling constants to two pions as the 
gi. It may be noted that the entire amplitude is proportional to s. The zero of the amplitude at s = is referred to 
as the Adler zero. Notice also that the higher terms involve the scalar masses and hence will vanish as the ^ oo. 
In the case of the linear-in-s current algebra term, the non zero result arose because the gi 's increase as mf . Taking 
the scalar masses to infinity is the same as integrating them out of the Lagrangian which results, as pointed out in 
the original paper [5l| by Gell-Mann and Levy, in a nonlinear sigma model. The magic cancellations in that case are 
very easy to see. Clearly they are more intricate in the present case. 

From the starting equation it is seen that the radius of convergence of the series in s is equal to the squared 
mass of the lightest scalar meson. To go beyond this point, in principle one should calculate all loop diagrams. A 
simple approximation is to identify the partial wave corresponding to the tree term with the K matrix amplitude. This 
gives results for amplitudes spanning a considerable range in s in reasonable agreement with present experimental 
indications. This was carried out for the SU(2) single M linear sigma model in 'u\ and for the SU(3) single M linear 
sigma model in |31| . 
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VII. FOUR POINT VERTICES 



We start by establishing the notations for the quadrihnear couphng constants involving the prediagonal fields: 

d^Vo 

Pi 

P2 



P3l = 



P32 



Pi 



d(l)\d(t)ld(t)\d(l)l 
dcj,'ld4>ld4>\d4>l 
d^'ld4>'ld<l>\d<j>l 
d<p'ld4>ld^'ld4>l 



= / ^ ^ \ (69) 

We find the following equations relating these quadrilinear coupling constants to the trilinear coupling constants in 
Eq. ([50l) by differentiating the second generating equation in (fT2|l three times with respect to pseudoscalar fields : 



P2 = Psi + — =- —=92 + —r2 + — <?3 + —rz 

a \/2a Vv3 



P32 = P-ii 



V\/3 



2 




2 






"71 


2 




2 






71 


2 


(f3 - f2)^ 








2 




2 






71' 


2 




2 


71^^ 




'71 



(3^ (3 f 1 2 1 2 \ V2 f 1 2 

1 / 1 



a 1/1 
P4 = -^P3i + ^ (^^.4 + -.4 + -.5 + -ra 

a f 1 2 1 2\V2/1 2, 

= ^^^^ - 7W [Ts'' ^ Ti' ^ Ts'' ^ TeV ^ T [Ts'' ^ 71'^ ' ^''^ 

Notice that the above equations were obtained by expressing five out of the six quantities in Eq. (j69p in terms 
of trilinear coupling constants as well as the sixth quadrilinear, pai. This shows that all the quadrilinear coupling 
constants cannot be obtained in terms of the trilinear ones. Nevertheless, as we will now see, the physical quadrilinear 
coupling constant, g can be completely expressed in terms of the bilinear coupling constants. Using the definition in 
Ea. (|62[) . we express the physical four point coupling constant in terms of the bare four point coupling constants as, 

/ d^Vo \ 

g = {Rt,)ai{Rt,)bi{.R-k)ci{R-k)di { a^^2^ a/j,l\ a^i.2\ a/ ) ' C^l) 

\ d{(l)i)Ad{(f>^)Bd{(l)i)cd{(l)^)D I 

which may be explicitly written as, 

g — cos'^OttPi ~ 4cos^6'^sin6'7rP2 + cos^0TrSin^6'7r(4p3i + 2^32) — 4cos0Tsin^07rP4 + sin^6'^p5. (72) 

Substituting Eg. ([70)1 into Eg. ([7^ and then using Egs. ([5^ and ([55)1 gives the formula for the quadrilinear coupling 
constant: 



1 



(a2 + /32)2 



r^{a^X^)n + 2a(3{X^),2 + /3'(^o')22) 



Noting a = {R-^)iiF^/2 and /? = {R-^)i2F^/2 we rewrite Eq.^ as, 



(73) 



9 = J2 (liR^')iD{XS)D.j{R.)ji + ^{R-')MX^)dj{R.)ji ) . (74) 
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In order to verify the cancellation of the s independent terms in Eq. ([64|) we should subtract half of Eq. ([74|) from the 
sum of the following two expressions: 

^ ':rE^i^-^)Al{XQ)AH{Lo)HG—^^ {R7r)ci{XQ)cK{Lo)KG 



GG 



71 



QP2 (^7r^)i£'(io)r>B(io ^)BA(^o)A//(-^o)i/G7:5T- — {Lq ^)gk{X^)kc{Lo)cf{Lp, ^)fj{Rtt) 

■7r^{R^^)lD{Lo)DE{Xl)EG—^. {XI)gf{Lq^)fj{R-!,)ji = T?2 (■^^^)l-D (^0 )£'7(^'r) Jl (75) 

>^^7r (^o)gG •^^tt 



r {R-r,)Al {XI)aH {Lo)hG (i?7r)ci {Xl)cK (^o) 



/fG 



[X^)dD (-^o)gG 



(^a^)GG 

= T^i^^^)^D{Lo)DE{Xl)EG—^ {Xa)GF{L'fi^)FjiR-^).n = -^-^{R^^)id{XI) d.j{Rii) .11 (76) 

It has thus been shown that the simple formula Eq. (p7)) holds near threshold in the case of massless pions for an 
arbitrary potential, Vq. 



VIII. SUMMARY AND DISCUSSION 



We have given a detailed treatment of a systematic approach to the study of a linear sigma model containing 
one chiral nonet transforming under U(l)yi as a quark-antiquark composite and another chiral nonet transforming 
as a diquark-anti diquark composite (or, equivalently from a symmetry point of view, as a two meson molecule). 
Some highlights of this work have been presented elsewhere [521 ■ The model provides an intuitive explanation of a 
current puzzle in low energy QCD: Recent work has suggested the existence of a lighter than 1 GeV nonet of scalar 
mesons which behave like four quark composites. On the other hand, the validity of a spontaneously broken chiral 
symmetric description would suggest that these states be (perhaps somewhat distorted) chiral partners of the light 
pseudoscalar mesons which are two quark composites. The model solves the problem by starting with the two chiral 
nonets mentioned and allowing them to mix with each other. Working with the SU(3) invariant version of the model 
it is seen that the four experimental inputs given in Eq. ((47| (note that the lighter 0~ nonet automatically has zero 
mass in the limit in which we are working) enforce a mixing whereby the light scalars have a large four quark content 
while the light pseudoscalars have a large two quark content. In addition, one light isosinglet scalar is exceptionally 
light (see Ea. (|ig)) ). 

Of the four experimental inputs just mentioned, there is a large uncertainty associated only with the mass of the 
"heavy" pion, the 7r(1300). It turns out that there is in fact some sensitivity to the precise choice of m[7r(1300)] so that 
this quantity is really being considered as a free parameter within the range of the quoted rather large experimental 
error. Thus the model parameters and predictions calculated in section IV are all displayed as functions of m[7r(1300)]. 
The effect of the not so large allowed variations in the mass of the ao(1450) are shown as error bars in these plots. 

In our treatment there are two parameters, associated with the masses and mixings of the SU(3) singlet pseu- 
doscalars, which describe the U{\)a anomaly in the effective Lagrangian. These parameters do not affect properties 
of the other particles and may be traded for the masses of the 7^(958) and one of the heavier candidates ry(1295), 
ry(1405), ?7(1475) or ?7(1760). The positivity of the eigenvalues allows only the last two candidates. For either of these 
it is noted in section IV that there are two solutions for the two quark vs four quark content of the r/(958). The 
presumably favored solution results in ry(958) with a mainly two quark content, while the less favored solution results 
in a mainly four quark content for the ?7(958). 

In sections V, VI and VII we gave a detailed proof that the low energy theorem for pion pion scattering holds in the 
present model with massless pions, for any choice of chiral invariant potential. The proof made use of the "generating 
equations" , stated in section II, to relate the four particle, three particle and two particle (ie mass term) vertices to 
each other. We carried out this somewhat lengthy calculation for two reasons. First, since the pion in the model has a 
non negligible, though small four quark content, one might wonder whether the theorem actually does hold. Second, 
it is expected to be useful to calculate the scattering amplitude, Eq. ([55|) in the resonance region, rather than close to 
threshold, as the theorem requires. 



19 



Clearly, there are a number of other interesting directions for further work. We plan to add mass terms in the 
same systematic scheme employed in section IV for selecting the most important chiral invariant terms. Mixing with 
glueball states and possibly other chiral nonets is also an intriguing possibility. Of course, an important ingredient 
to be taken into account would be the changes in the model parameters which result from unitarizing the tree level 
scattering amplitudes and comparing with the unitarizcd amplitudes with experiment. This was carried out for the 
2 flavor Cell Mann- Levy model in [6( and for the 3 flavor single M model in |3ll |. 
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APPENDIX A: SOME CORRECTIONS 

We have found the following minor corrections to ref [i^ : 

1. In Eq.(Al) the fifth term on the right hand side should properly read, (i2Tr(Af'M'^). 

2. In the sentence immediately following Eq.(Al), d2 should be added to the list of coefficients 
which are U{1)a invariant. 

3. In Eq.(19), the denominator of the argument of the "In" in the first term should read detM'l'. 

4. In the last line of Eq.(58) the left hand side should read (3^. 

5. In the last approximate equality in Eq.(60) the left hand side should read (3^. 



APPENDIX B: PARAMETER DETERMINATION 



Given the inputs: the pion decay constant, F^r; the mass of the ao(980), ttIq; the mass of the ao(1450), rua'] the 
mass of the 7r(1300), TOtt', the independent model parameters which don't involve the U(1)a violating terms can be 
successively determined (in the order given) by the equations: 



2d2 
4C2 

a 



mt, 



'"-a"'-a' 

■ml' - 

2 



ml, 

2 \2 



— ((rri^ — rn^,) 
\y a a I 



2d2 



-2(aeg) 
d2 

1 Fl 



4 1 + {P/af 



^[Ad2^{ml + ml,)f 
56(aeg)2 



1 

2^ 



C2 + 



8(aef)^ 



(Bl) 



The first equation tells us that d2 is positive for the experimental input masses. We take a and (3 to be positive. 
Then the fourth equation shows that 63 must be negative. Finally C2 and c% will be positive. 

Once the above parameters are determined, the parameters 71 and C3 of the U{\)a violating sector are obtained 
in terms of the mass of the ?7(958), m,,i and the mass of a suitable heavier 0~ isosinglet, m^2 as follows. First, 71 is 
found as a solution of the quadratic equation: 







= 572+T71 + C/, 
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U = r(l 



In addition, 



C3 



-36, 



j/o[m2 +m;2-yo(l + ^o) 



24yo 



(B2) 



(B3) 



Next we give the numerical values of the parameters for the central values of all the listed input masses except for 
TO[7r(1300)] which instead will take the typical value allowed by both the data and by the model, 1215 MeV. Table U 
shows the results for the parameters which arc not associated with the U{1)a violating part of the Lagrangian. 



C2(GeV2) 


9.64 xlO"^ 


d2(GeV2) 


6.32 xlO"^ 


eUGeV) 


-2.14 




42.1 


a(GeV) 


6.06 xlO-2 


/3(GeV) 


2.49 xlO"^ 



TABLE I; Calculated Lagrangian parameters:c2, d2, ef, c| and vacuum values: a, 13. 

Table HIl shows the calculated Lagrangian parameters associated with the U{1)a violating terms. Two "scenarios" 
associated with different identifications of the heavy which is the partner of the 7y(958) are shown (I assumes 7^(1475) 
to be chosen while II assumes 77(1760) to be chosen.) For each scenario, the two solutions (labeled 1 and 2) are shown. 





11 


12 


III 


112 


C3(GeV'') 


-2.42 X 10"* 


-2.42 X 10""* 


-3.44 X 10"" 


-3.44 X 10"" 


71 


5.38 xlO"^ 


2.53 X 10"^ 


8.69 xlO"^ 


-7.76 X 10-2 



TABLE II; Calculated parameters: cz and 71. 

Using these parameters we next list the mixing matrices for, respectively, the two 0" octet states, the two 0+ octet 
states and the two 0+ singlet states: 



0.925 0.380 
-0.380 0.925 



-0.496 0.869 
0.869 0.496 



0.711 0.703 
-0.703 0.711 



(B4) 



Similarly, the mixing matrices for the two solutions for scenario I of the singlet states are: 



I 1 : (^0"') = 



-0.671 0.742 
0.742 0.671 



, /2:(i?oi) = 



0.858 -0.514 
0.514 0.858 



(B5) 



Finally, the mixing matrices for the two solutions for scenario II of the singlet states are: 



-0.413 0.910 
0.910 0.413 



II 2 : (i?o ') 



0.974 -0.228 
0.228 0.074 



(B6) 
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